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Abstract. In this paper we determine orbital and linear stability of a class of spatially periodic 
wavetrain solutions with the mean zero property related to the Intermediate Long Wave equation. 
Our arguments follow the recent developments in m, da and [23] for the study of the stability of 
periodic traveling waves. 


1. Introduction 

One of the most fascinating phenomena given by nonlinear dispersive equations is the existence 
of solutions that maintain their shape and traveling with constant speed. Such solutions are caused 
by a perfect balance between the nonlinear and dispersive effects at the medium. In general, these 
solutions are called traveling waves and it is well known that the their existence has a very wide 
applications in fluid dynamics, nonlinear optics, hydrodynamic and many other helds (see pioneers 
works due to Boussinesq, Benjamin, Ono, Benjamin-Bona-Mahoney, Miura, Gardner, and Kruskal). 
Then, the study concerning the dynamics related to these solutions has became one of the important 
issues of the last decades for evolutive nonlinear partial differential equations. 

We can say that the initial impetus for the scientihc activity of these prohles was the inverse 
scattering theory (1ST) for the Korteweg-de Vries equation (KdV-equation henceforth) 

Ut “t“ Ux T (if )x ^xxx 0* 

One of the lessons learned by the 1ST is that the traveling wave with a solitary wave prohle, namely, 
u{x, t) = 'ip{x — ct) with c > 0 and 

lim i/>(0 = 0, 
kK+oo 

plays a central role in the long-time asymptotics of solutions to the initial-value problem associated to 
KdV-equation. Indeed, general classes of initial disturbances are known to solve into a hnite sequence 
of solitary waves followed by a dispersive tail. A companion result is that individual solitary waves 
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are orbitally stable solutions of the evolution equation. The exact theory of stability of solitary waves 
for the KdV-equation started by Benjamin in [15] (see also Bona IIZI) whose maturity was reached 
a decade ago with the works due to Albert jl], Albert and Bona [5|, Albert, Bona and Henry [6| and 
Weinstein |16|-|11|. Next, in papers due to Strauss at al. and Weinstein pa, pa, 051 were shown 
that not all solitary-wave solutions are stable. Both necessary and sufficient conditions for stability of 
the traveling waves solutions of a range of nonlinear dispersive evolution equations appear in various 
of the above references. 

In the last years, the study of stability of traveling waves of periodic type associated with nonlinear 
dispersive equations has increased significantly. A rich variety of new mathematical problems have 
emerged, as well as, the physical importance related to them. This subject is often studied in 
relation to the natural symmetries associated to the model (translation invariance and/or rotations 
invariance) and to perturbations of symmetric classes, e.g., the class of periodic functions with the 
same minimal period as the underlying wave. In the case of shallow-water wave models (or long 
internal waves in a density-stratified ocean, ion-acoustic waves in a plasma or acoustic waves on a 
crystal lattice), it is well known that a formal stability theory of periodic traveling wave has started 
with the pioneering work of Benjamin [151 regarding to the periodic steady solutions called cnoidal 
waves for the KdV equation. The waveform profiles were found first by Korteweg and de-Vries for 
KdV-equation. The cnoidal traveling wave solution, namely, u{x,t) = (pc{x — ct) has a profile given 
by 

</?c(0 = /^2 + (/^3 - /^2)cn^ ^ 

where cn{-;k) represents the Jacobi elliptic function called cnoidal associated with the elliptic mod¬ 
ulus k G (0,1) and are real constants satisfying the classical relations 

fdi < 1^2 < (dsi (di + 132 + /3s = 3c, ^ (1.2) 

P3 — Pi 


We recall that (pc satisfies the second order differential equation 


+ cMO - = Ac, ^ e 


(1.3) 


with = -I J2i<j f3if3j, and that the formula ( |l.l[ ) is deduced from the theory of elliptic integrals 
and elliptic functions. The existence of smooth solutions for (1.3) with a minimal period L, c G / C 
M —)■ (^c £ -^per([0)-b]) is determined from the implicit function theorem. The interval I in general 

depends of qualitative properties of (pc, for instance, for the property of mean zero, >c(0c^^ = 0, 
we have / = (0, J-oo) and for = 0 and <Pc(0 > 0 ^ G M, we have / = (^, J-oo). A first 

stability approach for the cnoidal wave profile (1.1) was began by Benjamin in [^ regarding the 
stability in L]) of the orbit 


= {^c{- + y) -yeR}, 


(1.4) 


by the periodic flow of the KdV equation. But only years later a complete study was carried out by 
Angulo, Bona and Scialom in [9] (see also 0 ). 

Recently, Angulo and Natali in [T3| (see also m have established a new approach for studying the 
stability of even and positive periodic traveling waves solutions associated to the general dispersive 
model 


Ut + 2uUx — (Mm)^; = 0, 


(1.5) 
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where M is a differential or pseudo-differential operator in the framework of periodic functions. M 
is defined as a Fourier multiplier operator by 

'Ji[g{n) = 6{n)'g{n), a G Z, (1.6) 

where the symbol 6 * of M is assumed to be a mensurable, locally bounded function on M, satisfying 
the condition 

< din) < 02(1 -I- (1-7) 

where mi < m 2 , \n\ ^ no, 0(n) > b for all n G Z, and a* > 0. One of the advantage of Angulo and 
Natali approach was the possibility of studying non-local evolution models in a periodic framework. 
For instance, let us consider the case of the Benjamin-Ono equation (henceforth BO-equation) 

Ui -\- uUx fff O 5 (7*S) 


with ffC denoting the periodic Hilbert transform and defined for L-periodic functions / as 


^fix) = ^p.v. 


rL/2 

r 7 r(x — 7)1 

/ cot 


/-L /2 

L L \ 


fiy)dy, 


(1.9) 


where p.v. represents the Cauchy principal value of the integral, we have that the Fourier transform 
of ‘Kf is given by the sequence {F£/(n)}„gz, where fff/(n) = —isgnin) fin). In other words, we have 
that M = 0-Cdx whose symbol is 6 *(n) = \n\. The periodic traveling waves uix, t) = ipcix — ct) for the 
BO-equation with minimal period L satisfies the following non-local pseudo-differential equation 

FC0C + C0c - R0C = 0) 


and they are given by 


dvr 


X = 


2 

senhij) 


L coshi'y) — cos(^) 


where 7 > 0 satisfies tanhi'y) = ^ (therefore the wave speed c must satisfy c > 27r/L). As an 
application of the theory in [13] , the authors obtained the first nonlinear stability result for the orbit 
generated by the wave ipc- 

In this paper, we are interested in studying the orbital and linear stability of a family periodic 
traveling waves for the physically relevant Intermediate Long Wave equation (ILW equation hence¬ 
forth) , 

Ut + 2uUx + S~^Ux - i75u)xx = 0, 5 > 0, (1-10) 

with u = uix, t) a L—periodic function and x, t G M. The linear operator 7s is defined by 

rL/2 

/ Ts,Lix - y)uiy)dy, 

'-L/2 


7suix) = ^P-v. 


where 


ryL(0 = -i^coth 

n^O 


(2^n.i/L 


\ L J 

Actually, the physical derivation of (1.10[) in a periodic setting requires that 


-L /2 


'-L/2 


uix)dx = 0 , 
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where we always can impose (1.10), because any non-zero mean could be removed by the Galilean 


transformation v{x,t) = u{x + 2'yt,t) — 7 , 7 G M. Hence, from the theory of elliptic functions (see 
Ablowitz, et al. [ 2 ]) we obtain that 

7iu{x) = —i^coth 




\L 


Moreover, for 6 ^ 00 , L hxed, we have (see 121 ) 


Jim 

S^oo \ Jj / 


which is the kernel of the Hilbert transform in (1.9). Therefore, the ILW equation (1.10) is the 


of the class of dispersive models (1.5) with exactly — t. 


natural periodic extension of the BO-equation (1.8). We note that the ILW equation is an example 


Now, one of our main objectives in this paper, it will be to hnd periodic solutions for (1.10) of the 


form u{x,t) = (pc(^ ~ ct) with the periodic prohle (pc having an mean zero and satisfying 

- C^Pc + Ac “ M<5 Ac = Ac 


rL 


( 1 . 11 ) 


where will be an integration constant given hy Ac = (Pc(x)dx. In section 3 we obtain, as 
a consequence of Theorem |3.1[ the following property associated to the pseudo-differential equation 


( 1 . 11 ): 


(FO) There is a smooth curve of even periodic solutions for (1.11) with the mean zero property, in 
the form 

c G / C M ^ Ac e i/”,([0, L]), n G N, 
all of them with the same minimal period L > 0 . 

By following the arguments due to Parker |32] (see also Nakamura and Matsuno in [H]), we obtain 


the following formula of even periodic solution for (1.11) with the mean zero property (see section 3 
below), 


cpc(x) := (Pc{L,d,k]x) = 


2K{k)i 


Z (^(. 


id); k 






( 1 . 12 ) 


where K{k) denotes the complete elliptic integral of the hrst kind, Z is the Jacobi Zeta Function and 
k G (0,1) (see notation section below). For hxed L and J, the wave-speed c and the elliptic modulus 
k must satisfy specihc restrictions. 

Other one focus of our study, it will be the dynamic of solutions of the ILW equation initially close 


to the mean-zero prohle (pc in (1.12), the stability of the prohle (pc- There are two common approaches 
to the stability question. Firstly, we can analyze the nonlinear initial-value problem governing the 
diherence between an arbitrary solution of the ILW equation and a given exact solution representing 
a wavetrain, the prohle ipc- In the hrst approximation, we assume that the diherence is small and 
we linearize the evolution equation. The resulting linear equation can be studied in an appropriate 
frame of reference by a spectral approach. To our knowledge, the linearized spectral approach has 
never been established for the ILW equation. A second approach to stability is the orbital stability, 
more exactly, we study the Lyapunov stability property of the orbit 


= {Vc(- + !/):!/ e 


( 1 . 13 ) 
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generated by the profile (pc- The study of the dynamic of the set consist in verifying that for 
any initial condition uo close to we have that the solution u{t) of ( 1 . 10 ) with m(0 ) = uq remains 
close to for all values of f G M. The specihc notion of “close” is based in terms of the following 
pseudo-metric dehned on a determined space W, namely, for f,gE W, 


d 2 {f,g) = inf 11 / - Trg\\w, 

r£K 


(1,14) 


with Trh{x) = h{x - 1 - r). The translation symmetry r enables us to form a quotient space, W/t, by 
identifying the translations r/ of each / G W. If we consider / and g as elements of IT/r, we obtain 
that ^2 represents a well-dehned metric on this set. Note that in W/t the difference u — (pc, between 
(Pc and the perturbed solution u, it will represent the most vital difference between two wave forms, 
namely, the shape. Again, according to our best knowledge, the orbital stability property associated 
to the prohle pc in (1.12) has never been established for the ILW equation in a periodic setting. 


Next, we shall give a brief explanation of our work. In fact, let us consider the new variable 

n(a;, t) = u{x -|- xt, t) — pdx), 


where u solves (1.10) and pc solves (1.11). Substituting this form in equation (1.10) and by using 


( 1 . 11 ) one hnds that v satishes the nonlinear equation 

2(vpc)x - cv^ - MsVx = 0 . 


Vt 


(1.15) 


As a leading approximation for small perturbation, we replace (1.15) by its linearization about pc, 
and hence obtain the linear equation 


Vt ^ dtr{Msv + cv - Ivtft). 


(1.16) 


Since pc depends only on x, the equation (1.16) admits treatment by separation of variables, which 
leads naturally to a spectral problem. Then, by seeking particular solutions of (1.16) of the form 
v{x,t) = e^^'ijj{x), where A G C, satishes the linear problem 

= A-/, (1.17) 

for £ := £c ,(5 denoting the self-adjoint operator 

£jc,s'■= + c — 2pc. (ITS) 


We recall that the complex growth rate A appears as (spectral) parameter. Equation (1.18) will only 
have a nonzero solution in a given Banach space Y for certain A G C. A necessary condition for 
the stability of pc is that there are not points A with Re(A) > 0 (which would imply the existence of 
a solution n of (1.16) that lies in E as a function of x and grows exponentially in time). If we denoted 
by a the spectrum of dx^, the later discussion suggests the utility of the following dehnition: 

Definition 1.1. (spectral stability and instability) A periodic traveling wave solution pc of the ILW 
equation (1.10) is said to be spectrally stable if a G iM. Otherwise (i.e., if a contains point with 
Re{X) > 0) Pc is spectrally unstable. 


We recall that as (1.16) is a real Hamiltonian equation, it forces certain elementary symmetries on 


the spectrum of a, more exactly, a is symmetric with respect to rehection in the real and imaginary 
axes. Therefore, it implies that exponentially growing perturbation are always paired with exponen¬ 


tially decaying ones. It is the reason by which was only required in Dehnition 1.1 that the spectral 
parameter A satishes that Re(A) > 0. 
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An similar spectral problem to (1.17) has been the focus of many research studies recently. For 


instance, if we restrict initially to traveling wave solution of solitary wave type, sufficient conditions in 
order to get the linear stability/instability has been established for many specific dispersive equations 
in Kapitula and Stefanov [36], in particular, the linear stability related to the generalized Korteweg-de 
Vries equation 

ut +{p + l)vFu^ + U:^xx = ^ p e N, (1.19) 

was obtained by using a Krein-Hamiltonian instability index to count the number of negative eigen¬ 
values with a positive real part. In the case of linear instability, Lin in |39| and Lopes in jlO] have 
presented sufficient conditions for general dispersive models. 

In a periodic framework, general spectral problem of the form 

JL'tp = X'tp 

has emerged, with J = and L a self-adjoint operator. Since J is not a one-to-one operator, 
classical linear stability results as in [28| can not be applied. To overcome this difficult, recently 
Deconinck and Kapitula in |21| (see also Haragus and Kapitula [20]) considered the similar problem 

( 1 . 20 ) 




in the closed subspace of mean zero. 


//„= /€T([0 ,L|); / f(x)dx = 0 


( 1 . 21 ) 

Thus, an specihc Krein-Hamiltonian index formula was deduced for concluding the linear stability 
of periodic prohle with a mean zero property. In particular, it was deduced the linear stability of 


periodic traveling waves of cnoidal type associated with the equation (1.19) for p = 2 (we also refer 
the reader to see Bronski, Johnson and Kapitula in 


and Deconinck and Nivala in [23]). We note, 
nevertheless, that for obtaining this specific result was necessary to know the periodic wave prohle 
as well as the knowledge of a specihc quantity of eigenvalues associated to the Lame problem 

—J- 6k‘^sn^{x; k)^ = 0$. 


Unfortunately, in our problem (1.17), this specihc type of information can not be established. 


We note that the spectral/orbital stability properties of periodic traveling waves in Hamiltonian 
equations that are hrst-order in time (e.g. the Korteweg-de Vries or the Schrodinger equations) have 
been very well-studied in recent years by using diherent approaches to those discussed above. See, 
for instance, Bronski and Johnson [TH], Bronski, Johnson and Kapitula [in]-[2n], Bronski, Johnson 
and Zumbrun [21], Deconinck and Kapitula [2S], Deconinck and Nivala [2E], Haragus and Kapitula 
[20] . Hur and Johnson [21], Jonhson [22]-[33] and Kapitula and Promislow [35] . 

In section 5 below, we use the approaches in Angulo and Natali [ID], Deconinck and Kapitula PI 
and Haragus and Kapitula [2D] for establishing the relevant result that the periodic prohle ipc in 
(1.12) for the ILW equation are linearly stable. By techniques reasons, we establish it result for c 
being strictly positive (see Remarks 4.1 and |5.2| below). 

Now, some informations for obtaining our linear stability result in section 5 for ipc in (1.12) can 
be used in order to conclude the orbital stability property of these periodic waves. Moreover, it 
property will be established for every admissible speed-wave c. Our approach, it will follow from a 
slight adaptation of the classical Lyapunov stability analysis established by Andrade and Pastor in 
[7]. In our case, the stability analysis will be based on the elliptic modulus k instead of c, such as 
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is standard in the classical literatnre, therefore we need establish a stability framework adapted to 
this new “speed-wave”. The energy space where the orbital stability property of the prohle (pc will 
be stndied, it is the following Hilbert-space, 


+00 


V0= {geLlML])-, 


W 


^ [l + 05 (m)]|^(m)p 


< CX) 


( 1 . 22 ) 


where ds indicates the symbol associated with M^. In section 6 , we briefly describe the main argn- 
ments for obtaining onr orbital resnlt of the prohle ipc by the periodic how of the ILW-eqnation. 

Onr paper is organized as follows. In section 2 we present notation and the dehnition of the Ja¬ 
cobi elliptic fnnctions. Section 3 is devoted to the existence of periodic waves having the mean zero 
property. In section 4, we present the reqnired spectral property associated with the linear operator 
(1.18) by following the arguments in na. In section 5, the linear stability of the periodic prohle (pc 
will be shown. To the end, in section 6 we establish our orbital stability result. 


2. Notation 


For k G (0,1), we dehne the normal elliptic integral of the first kind, 

X ^ 

/“ dt f dO ^. 

u(x-,k)= / — , = / , = Flip: k) 

J ^{l-t^){l-kH^) J Vl - k‘^sin^e 

with X = simp. The number k and (p are called the modulus and the argument, respectively. For 
X = 1 (</9 = |), the integral above is said to be complete. In this case, ones writes : 


K{k) = j 
0 


dt 

^{l-t^){l-~k^ 


dO 


\/l — k'^sin'^9 


Hence, K{0) = | and K{1) = J-cx). For k hxed, u = u{x;k) is a strictly increasing function of 
variable x (real). We dehne its inverse function by x = sn{u-, k) {snoidal function). Then, we obtain 
the basic Jacobian elliptic functions cnoidal and dnoidal, dehned by cn('u; k) = \/l — sn^(u-, k) and 
dn{u; k) = — k'^sn‘^{u; k) (see Byrd and Friedman [22] and Abramowitz and Segun |3|). Snoidal, 

cnoidal, and dnoidal have fundamental period 4:K{k), AKik) and 2K{k), respectively. Moreover, 
sn‘^{u; k) +cn'^{u; k) = 1 , kfsn^im, k) +dn^{u-, k) = 1 , sn(-u; 0) = sin{u), cn{u-, 0) = cos{u), sn{u; 1 ) = 
tanhipa) and cn(n; 1) = dn{u, 1) = sechipa). The Zeta Jacobi function, Zipa) = Zipa, k), it is dehned 
for M G M by 

U 

Z{u) = j 
0 

It is a function which is odd with fundamental period 2K{k). Moreover, Z{tt/ 2,k) = 0 and Z{mK) = 
0, para m = 0,1, 2,.... For u being a complex argument we refer the reader to formula 143.01 in [22] . 
In particular for n = ix, x G M we obtain 


dn‘^{x; k) 


m 

K{k) 


dx. 


Zipix, k) 


.snix] k') 

^ - 

cn(x; k') 


dn{x; k') 


iZ{x, k') 


TTX 

^2K{k)K{k')' 


with k' = Vl - 
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3. Existence of Periodic Waves. 

This section is devoted to establish the property (PO) defined in the introduction, more exactly, we 
construct a smooth curve of periodic waves with the mean zero property, c G I ^ ipc ^ ifpg^([0, L]) fl 
Hq, where the period L > 0 and the velocity c will have some specific restrictions. Our arguments 
will follow Hirota’s method, put forward in the works [H] and |12]. By convenience of the reader and 
from our stability approach to be established in sections 5 and 6 , we will review slightly the method. 

Indeed, let us assume the existence of / : C x M —)■ C, such that the profile 


, , . d 


In 


f{x + i6,t) 
f{x-i5, t) 


, (x, f) G M X 


it will satisfy equation (1.10), with being analytic in a specific rectangle R of the complex- 

plane. To simplify the notation, we define f+{x,t) = f{x + i6,t) and f-{x,t) = f{x — i6,t). So, by 
arguments in |32], there is a constant B, such that we have the bilinear equation 


iDf -\- B 


/+■/- = 0 , 


(3,1) 


with 


D"‘D'‘a(x, t) ■ b(x, t) := (9, - Si-)”(9i - t)b(x', 


(3,2) 


In addition, we can deduce from (3.1) that 

F(A,Px)/-/ = 0, 

where 

F{Dt,D^) = i ^Dt + smh{i6D^) + - B) cosh.{i5D^). 

Consider z = px + wt, where p, w G M will be determined later. Suppose that / has the following 
Jacobi Theta profile (see 0) 


f{x,t) = 93{z,q) := 1 + 2 


+0O 


cos( 2 nx) 


n=l 


+ 00 




for q = with r = i , where K'{k) = K[y/l — k'^) is the associated elliptic integral of the first 
kind. In general q = q{T) is the function called “nome” with Im(r) > 0. By substituting / at the 
identity (3.2), one has 

Fo93{2z, q^) + Fiq~^92{2z, q^) = 0 . 

Here, 02 represents the Jacobi Theta function of second kind. Moreover, one has 

+00 

Fm = F[2i{2n — m)w, 2i(2n — m = 0,1. 


In order to prove that f{x,t) = 93{z,q) is a periodic solution related to the equation (1.10), it is 
enough to prove that Fq = Fi = 0. To do so, it suffices to show that 


1 


- (w + f ) H'o - - AqB - 0 e T (w + ^) A[- ^A'[ - AiB - 0, 


6 


P 




P 


6 


p 




(3.3) 
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where 


Aq = Ao{p;q,S) 


+0O 

cosh(4nj95) = 6^{2ip6^ g^), 


n=—oo 


Ai = Ali(p;g,5) 


H-OO 

cosh[2(2n — l)p5] = q2d2{2ip5, q^) 


n=—oo 


and A[, i = 0,1, represent the derivative of the parameters Aq e Ai with respect to p, respectively. 
Next, we £x parameters p, q and 6 above. Solving the system in (3.3) we get 


B = B{p-q,6) 


52 yloA; - Al'oAll 


and 


w = w{p; q, 5) 


P p2 

(5 ^ (5 ■ AIoA; - Al'oAli 


p 

5 5 


dp 


{lii|W'(yl„,yli)]}, 


where W (Aq, Ai) = AqA'^—A'^Ai indicates the Wronskian of Aq and Ai. Now, if we use some standard 
identities concerning the Jacobi elliptic functions (see [3] and |22]), we deduce that f{x,t) = O^i^z^q) 
must satisfy the identity (3.2) provided that 


B = B{p;q,5) = -p^ 


9'l{2ip5,q) 
9i{2ip5, q) 


0'i(O,g). 


and 


w = w{p] g, J) 


+ 2ip^ ■ 
5 


9'^{2ip5, q) 
9i{2ip5, qy 


where 9i represents the Jacobi Theta function of hrst kind. 

Now, similar arguments can be used if one considers the slight change of variables 2 h-)■ 
case we see that 


-. In this 
2 


B = B{p;k,5) 


4 [9i{ip5,q{k)) 9[{0,q{k)) 


(3.4) 


and 


w = w{p; k, J) 



9[{ip5,q{k)) 

9i{ip5,q{k))' 


(3.5) 


where k G (0,1), k' 


_ _ 7vK{k') 

\/l — k'^, and q{k) = 
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Hence, we obtain that our hypothetic solution u becomes 




^3 ( -ipS),q{k) 


^3 ( ^{z+ ip6),q{k) 


(3.6) 


ip 

~2 


^3 ( \iz - ipS),q{k)\ 9'^ (^{z + ip6),q{k) 


^3 ( ^{z - ipS),q{k)j 03 \^^{z + ipS),q{k) 


it which represents a L-periodic function at the spatial variable with the natural choice of p = 2ti/L. 

Next, we obtain specihc restrictions on the parameter p, k and the minimal period L for u to be 
a smooth periodic function. Indeed, for k G (0,1) hxed, it is well known that the theta function 
6z{z,q{k)) has simple zeros at the points 


2; = (’11+2]^ + ^ '^1'^ 


So, the right-hand side of (3.6) possess inhnitely many isolated singularities which we need to avoid. 


To overcome this situation, it makes necessary to impose a convenient condition over the parameters 
p, 5 and fc, namely. 


0 < p(5 < —ilTT = 71 


K{k') 
K{k) ’ 


(3.7) 


(3,8) 


k' = i/l — k"^. To do so, it suffices to consider k G (0,1) satisfying 

25 K{k) ^ 

Our next step is to present a convenient formula for the solution u. Consider the parameters B 
and w satisfying condition in (3.4) and (3.5), respectively, then by using formula 16.43.3-[3] in (3.6) 
one has 


u{x,t) = 


iK{k)p 

71 

2K{k)i 


K{k) 


TT 


z — ipS); k 


( K{k) 


v ^ 


2 ; -|- ipS); k 


f2K{k), 

( —^—{x — ct — i5); k 


f2K{k) 


{x — ct + iS); k 


(3.9) 

(3.10) 


w 


para c :=-. Therefore, identity (3.9) determines a class of L—periodic functions which solves the 

P _ _ O 

ILW equation ( |1.10 ) with speed-wave c. Here, Z represents the periodic Jacobi Zeta Function (see 
section 2 above). 
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Next, we will determined an expression for c. Indeed, from the analysis above we obtain that 


w 1 . 9[{ip6,q{k)) 
p ~6 ei{ip6,q{k)) 


Thus, if we use formula 16.34.1 in [3], we get 



1 4:iK{k) 

6 ir~ 


^ V L ’ ^ fAi5K{k) 




sn 




-■,k 


(3.11) 


with sn, cn, dn denoting the Jacobi elliptic functions snoidal, cnoida and dnoidal, respectively (see 
section 2 above). Hence, for ^ = x — ct in (2.10) we obtain the periodic traveling wave solution ip^ in 


(1.12) for the ILW equation. Moreover, by construction one has that ipc ^ Hq. 


Next, by using formula 143.01 in |22], we can rewrite the prohle ipc iii function of the Jacobi elliptic 
functions snoidal, cnoida, and dnoidal: 


K{k) 


(2K(k)S \ 4fa 

L L '^) n Kit 

._(2K(k)x \ (2K(k)i \ (2K(k)i \ , / 

mk) ^j ■ j ■ =“ j ■ 


(3.12) 




1 — dn 


^f2K{k)x _^f2K{k)6 


■ 


k ■ sn 


V 




Figure 1 below, it shows the prohle pc with some specihc parameters of L,6 and k. 



Figure 1. Function pc iu (3.12) with L = tt, 6 = 1 and k = 0.5. 
































12 


STABILITY OF PERIODIC WAVES. 


Moreover, by using formulas 143.02, 161.01 and 120.02 in [22] at the identity (3.11) one arrives to 
the convenient formula for c = c{k), 




87r6K{k) 4:K{k) fA6K{k) 


h L^K{k') 




-■k' 


4:K{k) 

L 


cn 


nsKik) 


; fc' ■ dn 


f46K{k) 


k' 


sn 


r46K{k) 


■k' 


(3.13) 


Lastly, it follows immediate from condition (3.8) that for L and S fixed there is an interval (0, ki) C 
(0,1), with ki = ki(L, 6), such that v{L, 6,k) < 1 for all k G (0, fci). Therefore, we have the following 
existence result of periodic traveling wave for the ILW equation by depending of the elliptic modulus 

k. 


Theorem 3.1. For L and 6 fixed, there is ki G (0,1) such that for c 
have that the following smooth mapping 


c{k) defined in (3.13) we 


it is well defined. 

A = A(k) = i fj- 


k G (0, ki) —)■ fpcik) £ L]) F Hq n G N, 

Moreover, for every k G (0,/ci) we obtain that (pk = 
pl{x)dx. 


(3.14) 


Pc{k) satisfies (1.11) with 


In our analysis of linear and orbital stability of the prohle (pc{k) in sections 5 and 6 below, we need 
to determine the sign of the derivate ■^c{k). For arbitrary values of L and 5 this calculation becomes 
a challenge. By making many numerical simulations with fixed values of L and S we obtain that 
c = c(k) will always represent a strictly increasing function on the specihc interval (0, fci), and so 
we can assure the property (PO). For instance, the specific case of L = tt and h = 1 we obtain the 
following plots for the function c(k) and its derivate fi(k), respectively. 




Moreover, from the formula in (3.13) and some numerical simulations, we obtain immediately that 
ki in Theorem 3.1 has the approximation k^ ~ 0,944085037, and for k G (0,/ci) we have the basic 
condition in (3.8), v(7r, 1, k) < 1, and 


c(0) = lim c(k) 
k—>0+ 


-1.07462944, and 


lim c(k) = + 00 . 
k—kkr 
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We note that, there is a unique fco ~ 0.795178532 such that 

c{ko) = 0 and c{k) > 0, for all k G (3.15) 

therefore, the velocity c is negative on the interval (0, ko). 

The simulations for the cases L = 7i and differentes values of 6, showed a similar behavior of the 
functions c{k) and c'{k) as showed above. 


4. Spectral Analysis 


In this section, we start the analysis of the spectral problem (1.17) with L := Lc^s dehned in (1.18) 


The main idea for this study will be determine two specihc spectral properties for L, namely, that 
the kernel is one-dimensional with ker{L) = [;^</3c] and the existence of a unique negative eigenvalue 
which is simple. Since the operator L is non-local this analysis is not immediate. In this point we 
will apply the theory of Angulo and Natali put forward in [10] for studying the stability of periodic 
traveling waves for the nonlinear dispersive model (1.5). The initial obstacle for applying Angulo 


and Natali’s approach is that the periodic traveling wave prohle related to the equation (1.5) needs 
to be positive and satisfying the equation 

Mi/’ + c'lfj — = 0. 

Moreover, the wave speed ^ needs to satisfy > — infj.g]R 6*(r) in order to determine that M -|- <^ is 


a positive operator. In our analysis above (section 2), the traveling wave prohle of (fc in (1.12) has 
mean zero and the constant Ac in (1.11) is not zero. In order to overcome this difficulty, we shall use 


that the ILW equation is invariant by the Galilean transformation 


v{x, t) = u{x + 2'yt, t) — 7, 

for 7 being a real arbitrary value. The second obstacle is to determine the required spectral properties 
associated with the linearized operator L for arbitrary values of L and 6. So, by convenience in the 
exposition we shall restrict on a couple of specihc values for L and S, L = n and 5 = 1, respectively. 
However, numerical simulations enable us to conclude that for other arbitrary values of L and 6 our 
results remain valid. 

In what follows, we establish some preliminaries dehnitions and results due to Angulo and Natali’s 

in [13]. 


Definition 4.1. We say that a sequence a = (an)nez C M A m the class PF{2) discrete if 


i) an > 0, for all n G Z, 

ii) an^—m\(^n2—m2 ^n\—m2^n2—rai ^ 0? for Tli 71-2 and TTLi < 1/7.2, 

hi) m2 ^n\—m2^n2—mi ^ 0? if ^1 ^ 1^2; ^ 7772,772 ^ IRl, and n\ 7772- 

The dehnition above is a particular case of the continuous ones which appears in [1] (see also 
Karlin 133 ), namely, we say that a function : M —)■ M is in PF(2)-continuous if, 

i) g{x) > 0, for all x G M, 

ii) g{xi - yi)g{x 2 - 7/2) - g{xi - y 2 )g{,X 2 - yi) > 0, for xi < X 2 and 7/1 < 7/2, 

iii) strict inequality holds in (ii) whenever the intervals (xi,X2) and (7/1,7/2) intersect. 
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An sufficient condition for g belongs to PF(2)-continuous is for g to be logarithmically concave, 
namely, 

—log[g{x)\ <0, 

As an example of PF(2)-continuous functions, we have the prohle Qq{x) = sech^(a;), for p > 0, and 
for 0 < z/ < /i 

sinhiyx) 

Q{x) = —TT—r. 
sinn[jjix) 

Hence, the sequences (Qo(’^))nez and {Q{n))n^i belong to the class PP(2) discrete. 

The main theorem in uni is the following 


Theorem 4.1. Suppose that is an even positive solution of (1.11) with A = 0, namely, 

- 'i/'c = 0, 


such that ^ PP(2) discrete. Then the self-adjoint operator := M + ? — 2i/i^ possesses 

only one negative eigenvalue which is simple and zero is a simple eigenvalue with eigenfunction ■^4’,;. 
Moreover, its spectrum is bounded away from zero. 


the linear operator L in (1.18). 


Our focus in the following is to apply Theorem 4.1 in order to prove our main result associated to 


Theorem 4.2. Let L = tt and 5 = 1 and consider k G (0, fci), with ki defined by Theorem 3.1 
Then for ipc defined in (1.12) with c = c{k) 


we have that L in (|1.18|) is a self-adjoint operator 

d_ 
dx 


with a discrete spectrum and satisfying ker(£) = [^(Pc]- In addition, L possess a unigue negative 
eigenvalue which simple and the remainder of the spectrum is constituted by isolated real numbers 
which are bounded away from zero. 


Proof. Initially, from the specihc form of L we obtain from classical perturbation theory and spectral 
theory that £ is a self-adjoint operator with a discrete spectrum (see PI). 

Now, in order to simplify the notation, we denote 

N{k) := [ dx, R{k) := 

Jo 


N{k) 

L 


mi : = 


4K{k) f2K{k)6 f2K{k)6 


cn 




; • sn 


-; k'^ ■ dn ^ 




2 f2K{k)6 
m2 := sn I ———; 


4K(k) ^f2K(k)5 4671 Kik) , 2K(k) 

ms:= -1 - —and m, : = 


L2 K{k')' 


(4.1) 

(4.2) 

(4.3) 


In the following analysis we will leave the parameters L and 6 hxed, but arbitrary. Thus, from 


( |3.12[ ), ( |4.1[ ), ( |4.2[ ) and ( |4.3[ ) we get the expression 

dn^ (m4 • x; k) 


(Pc{x) = mi ■ 


1 — m 2 ■ dn (m4 ■ x; k) 


+ m3 


(4.4) 
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and, consequently, 


N{k) = 


dn"^ (m 4 • x; k) 


■ / - 2 -S 

io [l — m 2 ■ dn^ (m 4 ■ X] k)\‘ 


dx 


dn^ (m 4 ■ x]k) n 

2mim^ ■ I - -^^ 5 -^-^ dx + Lm^. 


/o 1 — m 2 ■ dn^ {mi ■ X] k) 

Next, by using formula 410.04 in [22] we deduce 


dn^ (m 4 • x] k) 
lo 1 — m 2 ■ dn^ (m 4 • x; k) 


dx = —■ dn 2 (C;fc) 


rui Jq 1 - m 2 ■ dn^ (C; fc) 

1 dn" (C; k) 


dC 


mi Jq 1 — m2 + m2 ■/c^ ■ sn^ (C; fc) 


dn^ (C; k) 


mi ■ (1 — m 2 ) 

2 

mi ■ (1 — m 2 ) 


dC 


dC 


0 l-a2-sn2(C;/c) 

TT ■ (/c^ — a^) ■ Ao(V’, k) 
2 \/ • (1 — q;2) • (q;2 — k "^) 


where 


a = 


m 2 




1 - m 2 


< 0 , m 2 7 ^ 1 , J) = sin 


-1 


a 


— k^ ’ 


(4.5) 


(4.6) 


(4.7) 


and Aq indicates the Lambda Henman function dehned by 


where 


Ao(V>, k) = -- [E{k) ■ F{ij, k') + K{k) ■ E{ij, k') - K{k) ■ E{J;, k ')], 

TT 


^(^) = / E{'ilj,k')=l \/l-{l-k‘^)siv?{e) dO 


and 


E{J;,k') = 


dO 


° W 1 — (1 — 


(4.8) 


(4.9) 


(4.10) 






















16 


STABILITY OF PERIODIC WAVES. 


Therefore, formula 410.08 in 123 enables us to conclude 


dn^ (m4 ■ X] k) 




'o [l — m 2 ■ dn^ (m 4 ■ x; k)~\' 


dx = — ■ 


[dn^ (C;fc)]' 


dC 


where 


and 


"^4 Jo [l — m2 ■ dn^ (C; A;)] 

1 r2K(k) [l-e.sn^{C,k)] 


Jo [l — m2 + m2 ■ • sn^ (C; k)~\ ^ 

2 /■A(fc) 1^1 _ /;;2 . gn2 


■ sn^ (C; k)\ 


dC 


dC 


U{a^,k) = 


V. = 


m4 ■ (1 - m2)2 Jq _ 

= - 

m4 ■ (1 — m2)^ 

+ 2-k‘^-{a^- k‘^) ■ n(a^ k) + (a^ - k‘^f ■ I/ 2 ] , 

k"^ ■ K{k) 71 ■ ■ Aq{'iP, k) 

^2 — 2\/a'^ ■ (1 — q; 2) ■ (a^ — 

[2-k^-a^-2-k^ + a^-{l- ■ K{k) 


(4.11) 


(4.12) 


+ a^-E{k) 
Statements ( 4.1[ )-( [4T^ , give us 


2 ■ (a^ — 1) ■ (fc^ — a^) 1 k"^ — 

TT ■ (2 ■ + 2 ■ — 3 ■ A;^) ■ k) 1 


(4.13) 


2'^a2 • (1 — Q;2) • (q;2 — A;2) 


/■ 


■ AV) 

+ 2-k'‘-(X - k^) ■ ^(a^ «:) + (a^ - ■ U] 


(4,14) 


2 ■ mi ■ m 3 

m4 ■ (1 — m 2 ) 


77 ■ {kJ — a^) ■ Ao(V’, A;) 
\/a2 ■ (1 — a^) ■ (a^ — A:^) 


+ Lmg. 


Next, by considering the specihc values of L = vr, 5 = 
0, 944085037), the existence of a = a{k) > 0 such that 

+ ca — i? = 0. 

In fact, one has 

— c + V + 4i? 


1, we obtain for all k G (0, fci) (A:i ~ 

(4.15) 

(4.16) 
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Moreover, by using that 


mm cpc[x = ^pc 

X&[0,L] 


we 


find via numerical simulations (see Figure]^ below) that 

a{k) > -ipc(k) (^^) , for all k € (0,ki). 


(4.17) 



Figure 2. Consider L = tt and 5 = 1. The continuous line gives us the behavior of the 
function a = a{k) for k G (0, /ci). On the other hand, the dashed line shows us the behavior 
of -Pc(k){L/2) for /c G (0,fci). 


Next, let us define ? = (,{k) by 


q ■.= c + 2a = a/c^ + 4i? > 0 


and we consider the translation function := 


By using (4.17), we conclude 0^ > 0. Moreover, 
periodic function one has that is also an even L—periodic function. Now, 


since (pc is an even L 
we claim that pc; satisfies equation (1.11) with A = 0 
/3 G M, it follows from (1.11) and (4.15) that 


Indeed, since 3v[s{(pc + P) = 'Ms(pc, for all 


-Mspc; - qpc; + Pi 


— M-Sp,; — Cpc; — 2apc; + PI 

-M5P>c - c{ipc + a) - 2a{ipc + a) + ((pc + af 
—M-sipc — cipc + (fl — [ca + a^) = 0. 


In what follows, we will verify that for all k G (0, fci), {pq{ji)}n^i G PF{2) discrete. We recall that 
such values of k satisfy the analytic condition (3.8). Applying formula 905.01 of [22] in (1.12), we 
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obtain 


2m 

(Pc[x) = — 


. ( 2m7r , , ( 2m7r , 

+00 sin ~ J +°° ^ 




4?! 


sinh 

+00 sinh 




miiKik') 
K{k) 

^2m7r(5^ 


, , m7iK(k')\ 


■*=1 sinh 


mnK{k') 


K{k) 


■ cos 


f 2rmix^ 


that is, 


^ +00 sinh 

, 471 

(x)=a+ — 2_^ - 

"1=1 sinh 


^2m7r5^ 


TmiK{k'] 

K{k) 


■ cos 


/ 2m'Kx\ 

v“^; 


So, the periodic Fourier transform related to the function 0^ is expressed by (/>^(0) = a and 

f 2mTi6\ 


b^[m) = — ■ 


sinh 


m7iK{k')^ ’ 
K{k) 


for all m G Z — {0}. 


Letting 


27i 6 nKik') 

h> := —— and /i : = 


K{k) ’ 


we obtain from (3.8) immediately that 0 < p < /i. On the other hand, by considering 

, sinh(z/x) 
smh(/za:j 

we see that 


^ [log(Q(a;))] <0, V x ^ 0. 


(4.18) 


(4.19) 


(4.20) 


(4.21) 


(4.22) 


Therefore, we obtain that Q G PF(2)-continuous (see |1]). In addition, we obtain the following 
specihc calculation to be used below, 

/ 2n5x\ 


hm^- 


L J 47rSK(k) 


sinh 


7rK{k')x\ L‘^K{k') ■ 

K{k) 


(4.23) 


Next, the following picture show us that the function 

271 


a{k) -— ■ t(L, 5, k) = a{k) — 2 ■ t(7i, 1, k) 

Jj 


for k G (0, fci), it is strictly positive mapping. 
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Therefore, we obtain for all k G (0, ki) the relation 

A'K6K{k) 271 X i\ 

Hence, the statements (4.21 )-( 4.24[ ) allow us to dehne a smooth function r : M —)■ M such that 

2nQ{x) 


(4.24) 


T(x : = 


V X G (— cxD, —1] U [1, +cxd) 


r(0) = a and r in (—1,1) such that r G PF{2) continuous. Therefore, we can conclude that 

G PF{2) discrete. 


Hence, from Theorem 4.1 we obtain that the linear operator — 20,^ admits exactly 

one negative eigenvalue which is simple and zero is also a simple eigenvalue whose correspondent 
eigenfunction is Lastly, we analyze the operator L. Indeed, since 

+ <^ — 20 ^ = M 5 + (c + 2 a) — {2ipc + 2a) = M 5 + c — 2 (pc = -L, 

then we obtain 


(4.25) 


ker(£) = ker(£,j_ 5 ) = 
This hnishes the Theorem. 


1 d M 


r d 1 


— 

idx J 


and ? 7 ,(£) = 1. 


(4.26) 


□ 


Remark 4.1. To study the behaviour of the function a = a{k) in (4.16) in order to determine that 
(4.17) holds for arbitrary values of L and S will induce enormous technical difficulties if we do not 
use numerical simulations for fixing values of L and S. Maple 16 enable us to conclude that (4.17) 
remains still valid for general values of L and 6 satisfying the analytic condition in (3.8). As a 
conseguence, the results in Theorem |4.2| can be established for general values of L and 6. 
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5. Linear Stability for the ILW-Equation 


In this section we establish our linear stability result for the mean zero traveling wave (pc in (1-12). 
For the convenience of the reader we will give some dehnitions and specihc sufficient conditions for 
obtaining our linear stability result (see | 23 ] and H). 

We start our study by establishing some dehnitions associated to the operator with £ = 

in dm nnd hI in dm. 


Definition 5.1. We define, 

(1) kr as the number of positive real eigenvalues (counting multiplicities) of the operator 

(2) kc indicates the number of complex-valued eigenvalues with a positive real part (counting 
multiplicities) of the operator dxJ^lno- 

(3) For B a linear operator with domain D{B), we define the linear operator Im{B)u = Im{Bu) 
for u e D{B). 


We note immediately from the later Dehnition, that since Im{L) = 0 then kc is an even integer. 
Next, for a self-adjoint operator A, we denote by n{{w,Aw)) the dimension of the maximal sub¬ 
space for which {w,Aw) < 0 (Morse index of A). Also, let A be an eigenvalue for dxJ^ and Ex its 
corresponding eigenspace. The eigenvalue is said to have negative Krein signature if 

fc-(A) :=n(('«;, (A:|^^)|g^w)) > 1, 


otherwise, if k~ = 0, then the eigenvalue is said to have a positive Krein signature. If A is a 
geometrically and algebraically simple eigenvalue for dxdi with eigenfunction fix then Ex = [fix \, and 
so 


0 , {fi\fi^\^^)fix) > 

1, (^A, < 0 


The total Krein signature is given by k^ := X]AeiK\{o} ('^)- Since Im(£) = 0 we obtain that k^ 
is an even integer. 


Definition 5.2. The Hamiltonian-Krein index associated to the operator dxkL is the following non¬ 
negative integer 

2CHam — kr kc A kj^ . 


Next, let us consider the quantity 

J=(£-^l,l). (5.27) 

We also note that for any / G ker{L)^ the quantity {L~^f,f) is always independent of h G J^~^f. 
Now, we denote by T> the 2 x 2—matrix given by 


D = 


(Aj iPc 'Pc) 

(T'Vc, 1) 


Then, from 


and 


(£-■ 1 , 1 ) 

we have the following results: 


Vc, 1 ) 

(T-M,!) 


(5.28) 
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Theorem 5.1. Suppose that ker(L) = Z/J 7 ^ 0 and D is non-singular we have for the 

eigenvalue problem in (1.20) the following relation 


3CHam = n{L) - n(J) - n(X)). 

We recall that ?7,(J) = 0 J > 0 and n(J) = 1 J < 0. An immediate conseqnence of Theorem 


5.1 is the following criterinm of linear stability. 


Corollary 5.1. Under the assumptions of Theorem 5.1, if kc = K = =0 then the periodic wave 

ipc is linearly stable. In addition, if Xnam = 1 then the refereed periodic wave is linearly unstable. 


Proof. The first part of the corollary is an immediate conseqnence of Theorem 2.7 in [21] (see also 
[in])- Now, if Xnam = 1 we dednce that k^ = I since kc and k~ are even nonnegative integers. Then, 
the spectral problem in (1.20) has a positive eigenvalne which able ns to dednce the linear instability 
of the periodic wave (pc. □ 


Next we establish onr linear stability resnlt associated to the periodic traveling wave (pc in (3.12). 


Since onr stndy will be based on Theorem 5.1, the valne of Xnam ninst be calcnlated. From Theorem 


4.2 we have that n(T) = 1. Next will prove that n(fD) = 1 and n(J) = 0 by considering the case of 
c being positive by technical reasons. For obtaining these qnantities we will need to calcnlate some 
expressions for J and det(D) in terms of the Jacobi elliptic fnnctions. More explicitly, we will obtain 
(see propositions below) the following explicit formnlas: 


J={£-^1,1) = 


d 

\ 1 

cL + 2 — 
dc 

/ p>l{x)dx 
Jo 


(5.29) 


and 






x)dx. 


(5.30) 


Thns, we will prove that ^ (p^(x)dx > 0 and therefore J > 0 and det(D) < 0. Hence, n(J) = 0 


and n(H) = 1. Therefore, from Theorem 4.2 and Theorem 5.1 we conclnde that Xnam = 0. Then, 


by Corollary 5.1 one has that the periodic wave (pc is linearly stable. Formally, we have the following 
linear stability result. 


Theorem 5.2. Consider c > 0. 
ILW eguation. 


The periodic traveling waver p>c in (3.12) is linearly stable for the 


The focus of the following propositions will be to show that J > 0 and det(!D) < 0. We recall that 
for convenience in the exposition we are considering L = tt and J = 1. We start by establishing the 
following main result. 

Proposition 5.1. For c > 0 one has ^ (pl{x)dx > 0. 


Proof. We start with the relation 


dc 



ipl{x)dx 


dk d 
dc dk 



dk 

dc 


Nfk). 


(5.31) 
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Thus, since c'{k) > 0, for all k G {ko,ki) (see (3.15)), we only need to establish the sign of N\k). 
Before that, it makes necessary to handle with the quantity N{k) in (4.5) for obtaining a convenient 
expression for our calculations. Indeed, from (4.18[) and Plancherel Theorem, we obtain 


+ 00 


N{k)= / ipl^^){x)dx = L ^ \(pc{m)[‘ 


Stt 


2 +°o 

E 

m=l 


sinh 


^2m7rS^ 


sinh 


rmiK{k' 

K{k) 


2 ’ 


(5.32) 


for all k G (0, ki). So, one can take the hrst derivative with respect to G (fco, ^i) in (5.32) to deduce 

( 

m ■ 


N\k) = 


IbvT' 


3 

E 


. _ f 2mTT5 \ 

2 



■ d 

'K{k')' 

-| > 

V L y 


r 1 K{k) )\ 


dk 

K{k) 



m=l 


sinh 


m7TK{k'] 


Since 


d 


K{k' 


K{k) 

[E{k) - K{k)] ■ K{k') + K{k) ■ E{k') 


n 3 


> . 


dk [K{k) _ 

we obtain immediately that 
This hnishes the proof. 


fc(fc2 - 1) ■iP(A;)2 


< 0, for all k G (0,1), 


N'{k) > 0, for all k G (fco, ki). 


(5.33) 

□ 


Remark 5.1. By using the proof of Proposition \5B\ and the numerical calculations made in Section 

(5.34) 


2 (see (3.15)) we see that N'{k) > 0 for every k G (0, ki) — {fco}- So, we have 


— ||(^c||^ > 0 , for every c 7 ^ 0 


Next we establish the formulas (5.29) and (5.30). 

Proposition 5.2. For every c > 0 we obtain J > 0. In particular, n(J) = 0. 
Proof. Since / = 1 G L]), for all s > 0, and M 5 (l) = 0 we get 

T(l) = M5(1) + c — 2ipc = c — 2ipc- 

Then, since ker(£) = T 1 and T pc, one has from (5.35) that 

1 = cL-^l - 2£“Vc- 

Thus 

c{T-H, 1) = (l,l)+2(£-Vc,l). 

Now, since c > 0, we get 

, 1 , L 2(£-Vc, 1) 

1 = - + ^ 

c c 


(5.35) 

(5.36) 


(6.37) 


Next, by differentiating identity (1.11) with regard to c we obtain 


d 


dc 


1 d 


^ IpFc = -Fc-t^\\p 


L dc 


(5.38) 
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Then, by applying the operator £ ^ at the equality (5.38) we deduce 


= Vc-“||y^ciP^ ^1- 


Hence, since (pc has the mean zero property we have 


d \ d 

—(fc, 1 ) = TT / ^c{x) dx = 0, 


dc 


dc 


and so, by combining (5.39) and (5.40) it follows that 


(5.39) 


(5.40) 


(5.41) 


Therefore, from (5.37) and (5.41) we arrive to the equality 

L/C CIC C 

Lastly, since > 0 (Proposition |5.l[ ), we get 

J=(£-il,l) =- 


d ,, 

cL + 2-y. 


(6.42) 


Thus, we obtain the formula in (5.29) and from the hypotheses on c and Proposition 5.1 we have 
immediately that J > 0. This hnishes the proof. □ 


Remark 5.2. From (5.42) we note that the requirement for c to be positive in Proposition 5.2 has 
only technical reasons. If we do not require c> 0, the study of 5 will depends on a “heavy” numerical 
calculations. Here, additional calculations in Maple 16 enable us to say that J > 0, for all c ^ 0 
{k e (0, fci) - {ko}). 

Proposition 5.3. For c > 0 we obtain det(2)) < 0. In particular, n{T)) = 1. 


Proof. We start by obtaining expressions for the elements of the matrix H in (5.28). Indeed, from 
([5)3^ and ([T42|, 




c(£-‘ll|,l) 




cL+2ly,\ 


(5.43) 


Hence, by using identities (5.36) and (5.43) and the fact that (pc G Hq we obtain 


(£ VcT'c) = 2^^ = 


2cL + 4— ll^^c 
dc 


(5.44) 
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Then, since J 7 ^ 0 (Proposition |5.2[ ) follows from (5.42), (5.43) and (5.44) that 


det(D) = - 


(T (po ^c) 


Vc 1)^ 


(£- 11 , 1 ) 






2cL + 4—||(pc||^ cL + 2 —||(y 9 c|p 

dc dc 


2 3dc 


(5.45) 


Therefore, we obtain the formula in (5.30) and from Propositions 5.1 5.2 we have det(D) < 0. This 
hnishes the proof of the Proposition. □ 


Remark 5.3. From (1.17) and the fact that L = we deduce that the positive and periodic wave 


(j)^ is also linearly stable. 

6 . Orbital Stability for the ILW-equation 

In the last section we have proved that the Krein-Hamiltonian index %na.m associated to the linear 
operator dxJ2 is zero, and thus the linear stability of the periodic traveling wave (fc was obtained. 
The next outcome of the theory is to obtain informations about the orbital stability of these periodic 
prohles. From the theories established in [28], [27], [2l] and [35l Chapter 5.2.2], we can deduce that 
ifc will be a local minimizer of a constrained energy, and so the orbital stability of these periodic 
waves is expected to be obtained provided we present a convenient global well-posedness result for 


the model ( 1 . 10 ). 


Now, the study of orbital stability can be based on an analysis of Lyapunov type (see 0 . P!-ra- 
IZ71-I2S1-I31-I1S!) and it will work very well when the integration constant Ac in ( | 1 . 11 [ ) is constant 
or zero. In the case of the integration constant Ac to be a function of the wave velocity c, as in our 
case, it does not seem to be immediate to apply this strategy. Thus, our following purpose will be to 
apply the recent development in Andrade and Pastor [7] to handle such situations and so to obtain 
the orbital stability of the prohle (pc for every c 7 ^ 0 (see Theorem 5.2 and Remark 5.2) 


We start our study by presenting the formal dehnition of orbital stability. 

Definition 6.1. We say that the periodic wave ipc in (1.12) is orbitally stable with respect to (1.10) 


in the space W in (1.22), if for all £ > 0, there exists 5 > 0 such that if uo G L]) fl W, 

s > 3/2, with 11^0 — Fc 
one has 


w < 5 and u{t) is the solution of ( 1 . 10 ) with m( 0 ) = uq, then for all t G 
inf \ \u{t) — ipc{- + s) 


|w < ^ 

Otherwise, the periodic wave ipc is said to be orbitally unstable. 


From Dehnition 6 T we have that some information about the global well-posedness problem for 
the ILW-equation need to be established. That is the focus of the following theorem. 

3 

Theorem 6.1. Consider uq G iLpg^([0, Lj). If s > -, then there is a unique u G C(M; i7pg^([0, L])), 
such that u solves the initial value problem 

Ut -|- 2uux — (M^m)^; = 0, {x, t) G M X M. 

m(0) = Uq. 


(6.46) 
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In addition, for all T > 0 the mapping data-solution 

«„ e //^,(|o,L]) ^ « e C([o,r]iff^,(|o,L])), 

it is continuous. 

Proof. See Abdelouhab et al. in [T]. 

The ILW equation has the following three basic conserved quantities, 


□ 


E_i{u) = u dx, Eo{u) = - 


and 


Ei{u) = - 


{'Msu)u dx — - 
o 


m dx 


m dx. 


(6.47) 


(6.48) 


Indeed, from Theorem |6.1 and density arguments we deduce that for all t, 

E_i{u{t)) = E_i{uo), Eo{u{t)) = Eo{uo), and Ei{u{t)) = Ei{uo). 

Moreover, the ILW equation admits the following Hamiltonian structure 

Ut = -2uUa, + {M.5u),c = dx{-u‘^ + M-su) = da;E[{u). 

Our purpose in the following is to describe Andrade and Pastor’s approach in the case of the 
ILW equation. We note from Theorem 3.1 that the wave-velocity, c, of our periodic waves in (1.12) 
may also depend smoothly on the elliptic modulus k, {k —)■ c{k), by equation (3.11)). Our stability 
analysis will be based on this parameter instead of the wave velocity parameter c, such as is standard 
in the classical literature. Therefore, we need to establish a stability framework based on this new 
“wave-velocity” parameter k. Thus, by following [7j and [28], we consider for every k G (0, fci) the 
following manifold in the space W, 


dc dA 

Sfc = <{ M e W; Mk{u) = Mk{(pk), where Mk{u) := —Eo{u) + —E_i{u) 


(6.49) 


where ipk = (Pc{k) and A{k) = ^ (pl(x)dx. We note that the strategy established in [^ is a 


generalization of the results in |3l] . The assumptions to obtain the orbital stability of in the sense 
of Dehnition 6.1| and by depending of the parameter k are the following: 


(FO) There is a smooth curve of periodic solutions for (1.11) in the form, 

fcGJ CM—iP” ([0, L]) fl Hq, n G N; 


(FI) /cer(£) = [£(pfc]; 

(F2) £ has an unique negative eigenvalue A, it which is simple; 
(F3) ^ ^ ^ 


L 


dk 


dk 


< 0 . 


Conditions (FO) — (Fl) — (F2) have been established for us in the Theorems 3.1 and 4.2 above. 
With regard to the condition (F 3 ), if we derivate the equation in (1.11) with regard to k is obtained 
the relation 

,'dipk\ dc dA 


dk' 


dk 


dk 


M'M. 
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Thus, by Proposition 5.1, Remark 4.1 and G Ho we obtain for every k such that c = c{k) ^ 0, 


L 


dipk\ f d^pk 


dk 


dk 


= ( 6 . 60 ) 


The main Theorem of this section is the following. 


Theorem 6.2. Let k G (0, ki) he fixed such that c = c{k) ^ 0. Then the periodic wave ipk = Tc{k) i 


m 


(1.12) is orbitally stable by the periodic flow of the equation (1.10) in the sense of Definition 6.1 


By convenience of the reader we give a sketch of the proof of Theorem |6.2[ The proof of the 
following two Lemmas follow from the ideas in 0 , 0 , 120 , and [34] . 

Lemma 6 . 1 . There is e > 0 and a C^—function, uj : Uflipk) i—>■ M, with 

Ue{(pk) ■= {u G W; \\u - ipkWw < £}, 

such that 

(m(-+ a;(M)), = 0, for all u^Uflipk)- 

Lemma 6 . 2 . ITe consider the conditions (Pq) ~ (-Pi) ~ (-P 2 ) — (-P 3 ) above, and the set 

Ak := {<!> G W; (<h, M'M) = ($, ^Tk) = 0}. 

Then, there exists a constant C > 0 such that 

(£$,<1)) >C'||<h||^, for all ^eAk. 


Now, for M G W we dehne the pseudo-metric 

p{u,ipk)-.= inf \\u- p>k{- + r)\\v;, 

r£[0,L] 

it which indicates the distance between u and the orbit generated by pk via the translation symmetry, 
namely, Qk = {Tk{- r) : r G [0, L]}. 

The following Lemma establishes the local minimal property of the prohle ipk on the manifold S*,. 


Lemma 6.3. We consider the conditions (Pq) — (Pi) — (P2) — (P3) above, and we define the functional 
ITk = El + cEq -|- RP_i. Then, there exist £ > 0 and a constant C{e) >0 satisfying 

Tfc(M) - Tfc((pfc) > C{e) ■ [p{u,ipk)]‘^, 

for all u G Uflpk) H S^.. 

Proof. Consider u G W. Since is invariant under translations one has 3^k{u) = 3^k{u{- + f)), for 
all r G M. Thus, it is sufficient to show that 

Tfc(M(- + a;(M))) - ITkiPk) > C ■ [p{u, <pfc)]^ 

that there is a constant Ci G M such that 

V := u{- + u{u)) -ipk = CiM'flipk) + 2/, (6.51) 

where y E T)k = [M'flpk]]'^ H [-^Tk]'^- Next, since Mk is also invariant under translations we can 
apply Taylor’s formula to obtain 

Mk{u) = Mk{u{- A u:{u))) = Mk{p>k) + {M'flipk),v) + 0(||t||^). 


where 00 is the smooth function obtained in Lemma [6.1[ Indeed, for u ^ Tik follows from Lemma 6.1 


(6.52) 
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Hence, since ?/ G one has {M'f,{ipk),v) = (M^((^fc), C'iM^((^fc)) = CiA^, where iV is a constant 
which is associated with the wave speed c. Then, since Mk{u) = Mk{(pk) we obtain immediately 
from (6.52) that 


C*! = O(ll^^llw)- 

Now, by applying a Taylor’s expansion to around u{- + oj{u)) = ipk + v we obtain 

Tfc(M) - Tfe(</9fc) = + o(||t||w)5 


(6.53) 


because of = 0 and 3^ki^k) = By using (6.51) and (6.53) we have {Lv,v) = {^Cy^y) + 

0(||t||^), and so we conclude 3^k{u) — S^ki^’k) = -^{^y^y) + o(i|T||w)- Next, since y G 'Bk, by Lemma 


6.2 there is (T > 0 such that {J^y,y) > C|||/||^. Thus, 

Jk{u)-3^kM>C\\y\\l + o{\\v\\l,), 


(6.54) 


where C > 0. Therefore, from (6.51) we deduce that for e: > 0, small enough, there is C* = C{e) > 0 
such that 


Jkiu) - Tfc(</3fc) > C||t||w > C[p{u, (pk)f 


This hnishes the proof. 


□ 


Proof of Theorem |6.2| , The proof of the result follows from Theorem |6.1[ Lemma |6.3| and a 
convenient adaptation of Theorem 3.5 in |28] (see also 0 ). By contradiction, we can select Wn ■ = 
Mn(-,0) G Ui{ipk) niL^g^([0,L]), s > |, and £ > 0, such that \\wn - 0, with 

supp{un{-,t),ipk) > e, 
t>o 


where Un{-,t) is the corresponding solution of (6.1). Let us consider £ > 0 satisfying Lemma 6.1 
From continuity of Unit) at f G M, we consider the smallest tn > 0 satisfying 


pi^Un(^'y , ipk') . 


(6.55) 


The following step in the analysis will be to determine the existence of > 0 such that «„«„(■, tn) G 
Sfc, for n large. This is exactly the point in the theory that we will apply the strategy in |^. Indeed, 
let us dehne /„ : M —)■ M, such that for n hxed, 

q ,2 f-L j-L 

fni^(^) -^k (Q^7tn,(' , tn)) ^ ' J (') tn^ | dx T O ■ ' J Un(^' ■, tn) dx . O Qn T Oihn- 

We note immediately that /n(0) = 0, > 0 and Mk{(pk) > 0. Thus, for all n G N there exists 

an > 0 such that fn{oin) = Mk{(pk)- In other words, there is (an)neN C M, satisfying 


Mk{anUn{-,tn)) = Mk{(pk), for all n G N, 


(6.56) 


that is, i^anUni^'^ tn))nGN C Sfc. 

dc dA 

Next, let 7k(u) := —EAu) and Bkiu) := —E^iiu). Then, since Eq and E^i are continuous 
dk dk 

mapping one has 7k{wn) —t 7k{<^k) =: S' 7^ 0, 7lk{wn) —^ Bk{pk) =■ h and Mk{wn) — Mk{<Pk), as 
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n —)■ + 00 . So, 

= lMk(anUn(-,tn)) - Mk(Wn)l = \Mk{Wn) - ^ 0 . 

as n —)■ + 00 . On the other hand, 

0 ^ \(y^7f^(^Wji^ ~|- 0iji'7.]^{Wri) (l7 S” ^)| — Qn S” l'^fc(R^n) 5^1 l^fc(R^n) ^1 

that is, 

Zn ■ 7 Oin'^kiO^n) t Q h. 


0, 


(6.57) 


Therefore, statement (6.57) gives us that (a„)„gN is a bounded sequence and therefore, modulo a 
subsequence, one has —)■ ao, as n —)■ +cxd. We will see that cio = 1. Indeed, from (6.57) we get 


(1 — ao) ■ [(1 + ao) ■ g + h] — 0. 


(6.58) 


Now, since 


h dk 

1 + - = 1 + = 1 + — 
9 ^k\9^k) 






= 1 + 


0 


7k i}Pk^ 


= 1 > 0 , 


we obtain that ao > 0. Therefore, since g follows from (6.58) that ao = 1. 

Next, we claim that 

p{Un{-,tn),anUn{-,tn)) -1 0, U -)■ +CX). (6.59) 

In fact, since p{un{-, tn), Pk) = there are G M and O 2 > 0 such that 

II ’I^n(') ^n) IIW ^ II (■) ^n) 7^fc(‘ T r'fi)||'W T ||7^fc(‘ T F,2)||'W ^ £ T ||7^fe(‘ T F}^)||'W 

that is, (||Mn(-,^n)||w)„gN i® ^ bounded sequence. Therefore, the convergence a„ —>■ 1 and the relation 

p(l^n(' 5 ^n) 5 Cy.n'^n^'i tn) ) ^ || ('; ^n) Cy.n'^n(^'i ^n)||w ^ |1 Cln| ' || ('; ^n) || Wj 

implies ( |6.59| ). Therefore, an application of the triangle inequality and (6.55) show that (a„Mn(-, tn))neN C 
Ue{pk)- Hence, from Lemma 6.3 we conclude immediately the convergence 


piotpUn^^' 1 tn}i Pk') 

Lastly, by using (6.59) and (|6.60|) we obtain. 


0, n —)■ +CXD. 


(6.60) 


„ Pip^n^'1 tn)y Pk') ^ pip^n^' 1 ^n); Cln^n(') ^n)) T pip^n'^n^': ^n); Pk) t 0, Tl ^ TOO, 


which gives us a contradiction. The proof of Theorem 6.2 is now completed. 


□ 


Remark 6.1. The positive and periodic wave (j)^ in (4.19) is orbitally stable by a direct application 
of the arguments in 0. 
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